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Abstract. The dynamics of an ensemble of bistable elements under the influence of noise and with global time-delayed 
coupling is studied numerically by using a Langevin description and analytically by using 1) a Gaussian approximation and 
2) a dichotomic model. We find that for a strong enough positive feedback the system undergoes a phase transition and adopts 
a non-zero stationary mean field. A variety of coexisting oscillatory mean field states are found for positive and negative 
couplings. The magnitude of the oscillatory states is maximal for a certain noise temperature, i.e., the system demonstrates 
the phenomenon of coherence resonance. While away form the transition points the system dynamics is well described by the 
Gaussian approximation, near the bifurcations it is more adequately described by the dichotomic model. 



INTRODUCTION 



GAUSSIAN APPROXIMATION 



Stochastic rate processes in bi- or multi-stable systems 
lead to many interesting phenomena observed in various 
scientific areas ranging from physics to social science, 
and have thus been studied for a long time. 

Here we consider a network of stochastically driven 
bistable elements whose distinct feature is a time- 
delayed coupling. The time delays are considered as uni- 
form and the network elements are assumed to be highly 
interconnected, so that the connectivity can be approxi- 
mated by a global all to all coupling. 

The dynamics of the network is numerically explored 
by using a Langevin model and analytically by using 
a Gaussian approximation and a dichotomous model, 
derived from the corresponding Fokker-Planck equations 
and Master equation, respectively. 



In order to theoretically study the dynamical properties 
of a globally coupled set of noisy bistable elements (with 
no time delay), Desai and Zwanzig 1 1] derived a hierar- 
chy of equations for the cumulant moments of the dis- 
tribution function from the multi-dimensional Fokker- 
Planck equation for the joint probability distribution for 
all elements. For large noise intensities, when the statis- 
tics of individual elements are approximately Gaussian, 
this hierarchy can be truncated. Applying this approach 
to our system yields the following set of equations for the 
mean field X and the variance M = N~ 1 — X) 2 , 



X = X -X 3 -3XM+eX(t-t), 
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(2) 



LANGEVIN MODEL 

The Langevin model consists of N equations, each de- 
scribing the overdamped motion of a particle in a bistable 
potential in the presence of noise and global coupling to 
a time-delayed mean field X(t ~t) =N~ l Jlf =1 x,(f - t), 



xi(t) = Xi (t) -xi(t) 3 +eX{t-x) + V2D£, (t), 



(1) 



where T is the time delay, e is the coupling strength of 
the feedback and D denotes the variance of the Gaussian 
fluctuations t, (t). 



DICHOTOMOUS MODEL 

To study the dynamics of a single bistable element with 
time-delayed feedback Tsimring and Pikovsky 1 2] used a 
dichotomous approximation. In the dichotomous approx- 
imation intra- well fluctuations of x,- are neglected, so that 
a bistable element can be replaced by a discrete two-state 
system. Applying this complementary approach to our 
system, allows us to express the mean field dynamics in 
terms of the hopping rates pn,2i which denote the prob- 
ability of a bistable element to change its state form — 1 
to + 1 and vice versa, respectively. The equation for the 
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FIGURE 1. Phase diagram of the Langevin model (crosses), 
the Gaussian approximation (dashed lines) and the dichoto- 
mous theory (solid lines). Double lines indicate hysteretic tran- 
sitions among phases. 
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FIGURE 2. The normalized coherence j8 of the oscillatory 
state at e = —0.2. as a function of the noise strength D, for 
networks of different size N as well as for the dichotomic 
theory. The coherence is given through j8 = Hco p /8co, where 
H is the hight of the dominant spectral peak at co p and 8co is 
the half-width of the peak. 



mean field reads 

X=pl2-P2l-(p2l+Pl2)X, (3) 

where the hopping probabilities are given by the Kramers 
transition rate |3], which in the limit of small noise and 
small coupling strength reads (cf. 01) 

_ v/2=F3eX(f-T) / lT4eX(f-r)\ 



PHASE DIAGRAM 

A numerical study of the Langevin model (Eq. [!) 
shows that the system undergoes ordering transitions and 
demonstrates multistability. That is, for a strong enough 
positive coupling the system exhibits a non-zero station- 
ary mean field and a variety of stable oscillatory states 
are accessible for positive and negative feedback. While 
the transition to the non-zero stationary mean field is sec- 
ond order (continuous), the type of the oscillatory transi- 
tions depends on the system parameters and can be first 
order (discontinuous), associated with hysteretic behav- 
ior, or second order. 

A linear stability analysis of Eq. yields the criti- 
cal coupling strengths £ and the frequencies CO of the ac- 
cessible oscillatory states. It is found that the oscillation 
periods can be arbitrary small meaning that a separation 
of time scales of the macroscopic mean field dynamics 
and the microscopic hopping dynamics of the constituent 
bistable elements is not possible. 

A comparison of the Langevin model with the Gaus- 
sian approximation and the dichotomous theory shows 
that while near the transition points the system dynam- 
ics is strongly non-Gaussian, it is in this regime well de- 
scribed by a two-state model (see Fig.Q, which allows 
for a complete analysis of the bifurcations of the trivial 
equilibrium. 



COHERENCE RESONANCE 

The considered system exhibits the phenomena of coher- 
ence resonance and array-enhanced resonance (see Fig. 
□ 

Both, Kramers random switching frequency p (see 01 
and the frequency of the oscillatory states CO, resulting 
from the coupling with the time-delayed mean field, de- 
pend on the noise strength, i.e., p = p(D) and CO = Co(D). 
Thus, the noise can be tuned so that the random hoping 
between the potential wells of the bistable oscillators is 
synchronized with the periodic modulation of the mean- 
field. This statistical synchronization takes place when 
co = np, where the coherence of the oscillatory states 
becomes maximal. 

For the N = I case the coherence resonance was ob- 
served by Tsimring and Pikovsky |2]. We observe that 
the resonance phenomenon not only persists in globally 
coupled networks with large N, but is enhanced, a prop- 
erty which was found in other systems and is sometimes 
referred to as array-enhanced resonance |4]. 



ACKNOWLEDGMENTS 

This work was supported by the Swiss National Science 
Foundation (D.H.) and by the U.S. Department of En- 
ergy, Office of Basic Energy Sciences under grant DE- 
FG-03-96ER14592 (L.T.). 



REFERENCES 

1. Desai, R. C, and Zwanzig, R., J. Stat. Phys., 19, 1 (1978). 

2. Tsimring, L. S., and Pikovsky, A., Phys. Rev. Lett., 87, 
250602 (2001). 

3. Kramers, H., Physica (Utrecht), 1, 284 (1940). 

4. Zhou, C, Kurths, J., and Hu, B., Phys. Rev. Lett., 87, 98101 
(2001). 



